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Abstract. We present the numerical evolution of a massive test scalar fields around a 
Schwarzschild space-time. We proceed by using hyperboloidal slices that approach future 
null infinity, which is the boundary of scalar fields, and also demand the slices to penetrate 
the event horizon of the black hole. This approach allows the scalar field to be accreted by 
the black hole and to escape toward future null infinity. We track the evolution of the energy 
density of the scalar field, which determines the rate at which the scalar field is being diluted. 
We find polynomial decay of the energy density of the scalar field, and use it to estimate the 
rate of dilution of the field in time. Our findings imply that the energy density of the scalar 
field decreases even five orders of magnitude in time scales smaller than a year. This implies 
that if a supermassive black hole is the Schwarzschild solution, then scalar field dark matter 
would be diluted extremely fast. 
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1 Introduction 

The scalar field dark matter (SFDM) is an alternative to ACDM that predicts flat density 
profiles in the center of galaxies and fits the abundance of substructure [1, 2], which are two 
of the problems nowadays of the cold dark matter model. The assumption of this model is 
that the dark matter is an ultralight scalar field particle, whose mass ms determines a cut off 
in the mass power spectrum of structures [1]. An interesting assumption in such analysis is 
that the scalar field potential is a cosh-like potential, that behaves as an exponential at early 
times and as a free field (quadratic potential) at late times. Moreover, it was found that the 
SFDM enjoys the same advantages at cosmic scale as the standard lambda cold dark matter 
model [1]. 

Initially, within the SFDM model, fully general relativistic stationary solutions were 
proposed to explain the flatness of rotation curves [3]. Later on, the assumptions relaxed to 
the newtonian limit of such general relativistic models. By the time, the fluid dark matter 
made of scalar fields was proposed as an alternative galactic dark matter model [4]. On the 
other hand, the assumption of time independence was also relaxed and scalar field dark matter 
halos were proposed to be gigantic Oscillatons, that is, time dependent fully relativistic scalar 
field solutions to the Einstein-Klein-Gordon system of equations [5, 6]. On the other hand in 
order to study late time stages of the structures, the model relaxed to the Newtonian low- 
energy limit ruled by the Schrodinger-Poisson (SP) system of equations, where the turnaround 
point and speed of a structure formation was studied [7]; the condition to achieve quick galaxy 
formation is that the mass of the boson is again ultralight {uib ~ 10~^^eV) which corresponds 
well to the cut-off needed in the mass power spectrum [1]. In [8] it was also shown that the 
scalar field gravitational collapse tolerates the introduction of a self-interaction term in the 
potential, which makes the model to seem quite like a self-gravitating Bose-Condensate whose 
equilibrium solutions are attractors in time. In [9] the non-spherical collapse was studied and 
shown that again spherical equilibrium solutions are attractor profiles and in [10] it was 
shown the head-on collision of two structures and interference patterns were found during 
the collision. 

In this paper we explore the evolution in time of a test scalar field on top of the space- 
time background of a Schwarzschild black hole, and analyze the time development of the scalar 
field amplitude and its energy density in order to impose restrictions on the scalar field dark 
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matter model or the system Black Hole plus scalar field dark matter. An important property 
of the supermassive black hole plus dark matter system is the asymptotic nature of the space- 
time. As far as we can tell, black hole candidates are considered to be asymptotically flat, 
this is the reason why here we consider a Schwarzschild black hole as the model of black hole. 
Sporadic models like that in [11] assume a non- asymptotically flat space-time modeling the 
background of a galaxy, which could eventually explain high velocities of test particles in the 
center of a galaxy and at the same time the flat rotation curves. This sort of models assume 
there is a non-zero background energy density in the whole space-time which is a property 
inconsistent with asymptotic flatness. 

The standard approach appropriate for the study of the common dark matter candidates 
of the wimp type, modeled by collision-less matter near supermassive black hole's region is 
focused on two major streams: i) explore the possibility that current supermassive black 
holes are the result of the growth through accretion by seed black holes of intermediate 
mass of about 10^ to 10^ of baryonic and dark matter [12], and ii) the contribution of 
dark matter to the growth of supermassive black holes. Concerning the second stream, some 
results indicate that only about 10 % of their mass is due to the accretion of dark matter (e.g. 
[13]), and more mature models involving the study of the Newtonian phase space [14] indicate 
that the time-scale for the accretion of collision-less dark matter to contribute significantly 
to a super massive black hole mass is too long (see e.g. [15]). Also, important bounds on the 
collision-less dark matter central density have been proposed that allow the observed masses 
of central black holes without runaway instability of black holes, based on the study of orbits 
traveling on unbounded orbits with respect to the central object [16]. More recent results 
show the importance of pressure of dark matter in the accretion rates, where it was shown 
that pressure itself is able to freeze the black hole's growth through dark matter accretion 
[17]. 

In this paper we proceed in a similar way, that is we study the evolution of a scalar 
field acting as dark matter around a supermassive black hole and study its behavior. We also 
make use of the approximation that considers that the dark matter is diluted and thus its 
evolution affects the black hole space-time in a minor way, thus allowing one to consider the 
dark matter to be a test matter field as assumed in those of the aforementioned studies. To 
our purpose it is particularly important noticing that a property of a test massive or massless 
scalar field evolving on a Schwarzschild space-time, is that it can only move toward two ends: 
the inside part of the black hole and future null infinity (scri +). That is, if the scalar field 
evolves around black holes, either supermassive black holes, or less massive black holes, we 
wonder how the scalar field behaves in time, whether it remains near the black hole, or is 
accreted by the black hole or it leaks toward scri +. Previous studies invoke the quantum 
nature of the scalar field to understand the accretion rate of SFDM into a black hole, that 
is, given that the scalar field is ultralight, it has a Compton wave-length much larger than 
the Schwarzschild radius and therefore the black hole will have a restricted cross section for 
such a scalar field (see for instance [18]). Unlike such approach, in this paper we consider the 
scalar field to be classical as assumed in SFDM models at cosmic scales, from an effective 
Lagrangian: 

C = -R+{V^f + V{^), (1.1) 

where R is the Ricci scalar of the space-time, ^> the scalar field and V{^) its potential. The 
variation of such Lagrangian with respect to $ reduces to the Klein-Gordon (KG) equation, 
that rules the evolution of the scalar field, and since we work on a fixed background space- 
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time there is no need to variate the Lagrangian with respect to the metric, which would 
imply Einstein's equations. 

Thus, we study numerically the behavior of a scalar field playing the role of SFDM on 
top of a black hole space-time. In order to study the evolution of the scalar field we assume 
the following conditions: i) we set up a fixed Schwarzschild black hole space-time using 
hyperboloidal slices, and compactify the slices using coordinates such that our numerical 
domain contains scri +, ii) the compactification of the spatial coordinate implies a singularity 
of the metric at future null infinity, thus we define a conformal space-time and solve the KG 
equation on top of such conformal space-time, iii) we consider the scalar field to be a test field, 
that is, we do not solve for the evolution of the geometry of the space-time; this assumption 
works as a first approximation for a very diluted matter field, iv) we solve the Klein-Gordon 
equation as a hyperbolic initial value problem, v) study the late time behavior of the scalar 
field and its independence on the initial data, so that we can determine an accretion-leak rate 
of the scalar field, vi) we trigger the evolution of the scalar field by setting up initial data 
consisting on a time-symmetric shell of scalar field that presents both, an ingoing and an 
outgoing pulse, vii) based on the study of the evolution in time of the scalar field amplitude 
we build up our conclusions. 

Due to the considerable amount of calculations involved in the method we use, the units 
we assume in the paper are such that G = c = 1 and the mass of the black hole M = 1. 

The paper is organized as follows, in section 2 we describe the construction of coordinates 
describing a Schwarzschild black hole using constant mean curvature slices and scri-fixing 
compactification, in section 3 we describe the method used to solve the Klein-Gordon on the 
conformal space-time. In section 4 we present the evolution in time of the scalar field and 
interpret our results in terms of SFDM. Finally in section 5 we present some conclusions. 

2 Description of the compactified Schwarzschild space-time 

An appropriate coordinate system of the Schwarzschild space-time to study the propagation 
of scalar fields can be achieved using hyperboloidal slices, because such slices reach future 
null infinity instead of spatial infinity, which is a natural boundary for a wave-like process, 
including electromagnetic fields and gravitational radiation. Additionally we require the 
slices to penetrate the event horizon of the black hole in order to track the accretion of the 
scalar field up to the interior of the black hole. 

On the other hand, it is possible to compactify the space-time in order to work on a 
numerical domain that contains scri +, for which it is necessary to compactify the spatial 
coordinates in an appropriate way such that the space-time is regular there, at future null 
infinity. Usually, when a spatial coordinate is being compactified, the metric is singular at 
the new spatial infinity, and it also happens in the present case when hyperboloidal slices are 
compactified. Among others [19, 20], a known way of fixing this problem is the definition of a 
conformal metric that uses an ad hoc conformal factor that regularizes the singular terms in 
the metric resulting from the compactification [20, 21]. The idea behind scri-fixing conformal 
compactification for spherically symmetric space-times is the following: i) use hyperboloidal 
foliations whose slices reach scri +, ii) compactify the radial coordinate, iii) regularize the 
metric with an appropriate conformal rescaling. 

Assuming a spherically symmetric space-time is labeled using coordinates {t,r,6,ip), 
a well known process to foliate the space-time using hyperboloidal slices starts by defining 
a new time coordinate t = t — h{f), where h{f) is a function called height function. This 
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transformation has the advantage that keeps the time direction invariant. A compactifying 
coordinate is introduced in the form r = fi, where Q = Q{r) is a function that compactifies 
r and at the same time is a non-negative conformal factor that vanishes at scri +; in the 
present case we choose Q = 1 — r and then scri + is located at r = 1. We rescale the space- 
time using a conformal transformation g = Q'^g so that the conformal metric g is regular 
everywhere because the conformal rescaling compensates the singularities in the metric due 
to the compactification of the radial coordinate [21]. 

For the construction of h{f), constant mean curvature slices (CMC) are used [22]. The 
mean extrinsic curvature k of the initial slice t = is given by k = V^n^ = --^=dfj,{y/—gn^), 

where is a time-like unit vector normal to the spatial hypersurface and positive k means 
expansion. In this case, we adopt the requirement that the mean extrinsic curvature k is 
constant. This condition allows the integration of the equation above. For the Schwarzschild 
metric 

ds'^ = -a^{f)dP + a2(f)df2 + f2^(i02 ^ Q^^2^^ (2.1) 

where a = 1/a = ^Jl — it is found that, after solving for h' the condition k = constant, 
the final version of the conformal metric is [22]: 



= _ _ ^^ O^dt^ _ ^MZL^l^dtdr + ^dr2 + r\de^ + sin^ dd^''), (2.2) 
\ r J "v) F'^{r) 

where 



P{r) = n^P{r) = 

The values of k and the integration constant C are restricted such that -P(r) is real; the units 
of C are those of M^. The values of these constants used to obtain our are k = 0.4 and 
C = 2 in all our calculations, and it has been shown that for the massless case the results 
are independent of the values of these constants [21]. 
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3 Klein-Gordon equation: solution and numerical aspects 

The Klein-Gordon equation for a scalar field 6^ reads: = 0, where 0(hx = 

a(pT 

-^^df^[^—ggf^'^d,y(f>T], and we consider the potential to be of the form V = ^m^|i^r|^+f [f^Tl^, 
where nis has units of mass. The KG equation in the conformal metric reads 



n\(j)T - -R4'T 

6 



(m|0r + A0|) = 



6 ^ 



(m 



2 O- 



1)T + 



a4) 



0, (3.1) 



provided the relationship between the physical scalar field (pT and the conformal scalar field 



is (^T = 4>t/^- Here R 



120 



(r -|- (2r — 1)) is the Ricci scalar of the conformal metric 



and □ = V^V^ corresponds to the conformal metric. The case for which equation (3.1) is 
conformally invariant corresponds to the zero mass case niB = 0, however, as we want to 
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study the behavior of the test scalar field on a fixed Schwarzschild background, we solve 
the conformal KG equation that contains the mass parameter diff^erent from zero and avoid 
the future null infinity boundary at r = 1 in our calculations and we only approach such 
boundary in order to consider a big spatial physical domain. 

The scalar field is separated in the form (prit, r, 9, ip) = (^(i, r)Yim{0, ip) for a given value 
of I, in order to explore the behavior of angular components of the scalar field in the process. 
As mentioned above, since the conformal KG equation (3.1) is singular at r = 1, we do not 
have the chance to solve the equation in the whole spatial domain and we are forced to set 
up an artificial boundary close to scri + approximately at a distance 1000 in units of the 
mass of the black hole. We also implement an expansion of the factor in (3.1) in a power 
series Xl^o such boundary we apply outgoing wave boundary conditions, and since 

the boundary is pretty close to scri + the ingoing velocity is nearly zero and no numerical 
incoming signals affecting our calculations were observed. 

We implement another boundary inside the black hole's horizon at Vmim a method called 
excision. Notice first that the black hole's horizon according to (2.2) is located at r = 2/3. 
For our calculations we impose the boundary at rmin = 0.6, where we know the light cones 
point toward the singularity and no signals are propagated outward, which also involves 
spurious signals due to numerical errors, this is why there is no need to impose boundary 
conditions at this excision boundary. 

We solve (3.1) as an initial value problem using a first order formulation and usual ADM 
gauge functions. In terms of the line element (2.2) we can read off the gauge in terms of the 
ADM-like metric df = {-a^ + j^^'^)dt^ + + ^"^dr"^ + r'^{de'^ + sin^ Odip^) and obtain 

the following gauge and metric functions 



a^{r)h'{r){n-rn') 



which are the final metric functions we use to solve the KG equation. For numerical con- 
venience we solve the KG equation as a first order system, for which we define first order 
variables vr = ^dt(j)T — ^f3dr(j)T and ip = dr<j)T- In terms of these new variables the KG equa- 
tions splits into the following system of first order equations (see the appendix for details) 



dttp = dr ( ^TT + Pip 



dt4>T = "tt ^tp, 
7 

drCpT = (3.3) 

where the first two equations correspond to the evolution of the first order variables, the third 
one is the definition of vr that we use to reconstruct the value of (pT, and the last equation is 
a constraint that has to be satisfied during the evolution. We have included the contribution 
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of the angular part of the scalar field. System (3.3) is a symmetric hyperbolic system, which 
implies that we start up with a well posed initial value problem [23]. 

We integrate system (3.3) numerically using a finite differences approximation on a 
uniformly discretized domain, we use the method of lines with sixth order accurate stencils 
along the spatial direction and a fourth order Runge-Kutta time integrator. In all our runs 
we use a constant Courant-Priedrichs-Lewy factor (CFL), such that At = CFL x Ax, where 
At and Ax are respectively the time and spatial resolutions of our discretized domain. In 
order to set initial data consistent with (3.3) we restrict to a particular set of cases. That 
is, for the case A = eq. (3.3) is the same as for (p{t,r) and I = 0,1,2,..., because the 
spherical harmonics are factored out, for this reason initial data for (f){t, r) and any value of 
/ is consistent. On the other hand, for the case A 7^ 0, eq. (3.3) depends on the spherical 
harmonics due to the (jij. term; in order to simplify our parameter space we choose to provide 
only initial data for ^(t, r) and only the case I = for which spherical harmonics are trivial 
and the initial data are consistent. We use initial data corresponding to a time-symmetric 
Gaussian profile for the scalar field because it includes both, a mode moving toward the hole 
and an outgoing pulse of scalar field. The initial data used then read 

</.(0,r) = Ae-^^-^")'/'^', 
^(0,r) = -2i!l^<^(0,r), 

7r(0,r)=0, (3.4) 

and explore the parameter space for the amplitude and width of the Gaussian profile, while 
we keep tq fixed at rg = 0.8 which in the physical radial coordinate means fo = 4, and is 
well out of the event horizon. The value of a is extended to cover various length scales, given 
that for r > 0.8 the physical coordinate f grows rapidly as f = 

4 Results 

In this part, unless otherwise stated we use initial scalar field profiles with m?^ = 0.0, 0.1, 0.2. 
In order to explore the parameter space, we use various values of the amplitude A = 
0.01,0.1, 1.0 and different widths of the Gaussian pulse, which in physical units correspond 
to o"! = 0.5, (72 = 1 and <T3 = 5 (in units of the mass of the black hole) to the right from 
ro = 0.8. This range of parameters allows us to cover length ranges that involve potential 
Compton wave-length related effects of interaction and reflection suggested recently [24] , that 
may involve reflection and absorption effects. 



The evolution is tracked during a time window in which our results converge and that 
also allow to track the amplitude of the scalar field from the quasinormal mode ringing stage 
(in the massless scalar field case these are actually the quasinormal mode oscillations of the 
black hole) where the amplitude measured at various distances from the horizon behaves as 
~ Qi(i^R+i(^i)t ^ polynomial tail decay stage of the amplitude, during which the time 

dependence of the amplitude is of the type \4>\ ~ with p < 0. 

We measure the amplitude of the scalar field at various distances from the black hole's 
horizon. In Fig. 1 we show the evolution of its amplitude, and find as expected an oscillatory 
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Figure 1. (Left) We show the amphtude of the massive scalar field in time measured at various 
distances from the black hole's horizon for I ~ 0, 1,2 and = 0.1. The time dependence shows 
an oscillatory stage similar to that found for the massless case followed by a polynomial-like decay 
rate ~ with p < 0; each line indicates the scalar field amplitude measured at a different distance 
from the black hole in units of the mass of the black hole. (Right) We show the exponent p for the 
various cases. We show that the decay rate of the scalar field is bigger with bigger values of /, that is, 
the amplitude decreases faster during the tail-like stage when non-trivial angular components of the 
scalar field get involved. If this scalar field is to be considered to be the dark matter, this indicates 
that the amplitude decays faster if the scalar field has angular components. 



lapse followed by an approximately polynomial tail decay (the exponent does not stabilize 
around a constant value) for various locations of the detectors and three different values of 
I. We also show the exponent of the polynomial decay p in each case. We notice that the 
amplitude decays with decreasing exponents between —6.5 < p < —4, depending on the 
location of the detector for the / = case. For the cases with non trivial angular component 
the decay is even faster, for instance for the / = 2 case the exponent lies between —9 < p < —5 
even if such exponent has not yet stabilized. In these results we only show the time window 
within which our calculations converge at least with fourth order, which is the accuracy order 
of the time integrator as shown below. 
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Figure 2. In these figures, we sliow tlie amplitude of the massive scalar field measured at 1000 from 
the black hole's horizon, and for different values of the amplitude A and widths a of the Gaussian. We 
show the evolution for three values of the angular term I = 0,1,2 (from top to bottom respectively) 
and scalar mass rn^ = 0.1. We found that the tail behavior of the amplitude is about the same. This 
implies that scalar field decays at same rate independently of the profile parameters. 



In Fig. 2 we show the time dependence of the scalar field amplitude for various values of 
the parameters of the initial data. The finding is that the tail behavior is independent of the 
width and amplitude of the initial gaussian we explored, which indicates a universal behavior 
which does not seem to depend on the initial data. That at the tail stage of the amplitude 
the exponent of the decay is independent of the amplitude and width of the Gaussian profile 
allows us to extract conclusions for an initial SFDM packet represented by the scalar field. 

In Fig. 3 we show the amplitude of the scalar field for different values of the boson mass 
m^, measured by a detector located at r = 1000 for / = 0. We were unable to explore the 
parameter space for bigger values of because the tail-like decay appears at considerably 
small amplitudes near the round-off error threshold of our calculations, and we were unable 
to achieve convergence at such regimes, however this figure suggests a trend: the bigger the 
values of the faster the amplitude of the scalar field decays; this is consistent with the 
fact that the massive term in our evolution equations is always switched on and acts as a 
dissipative term. 
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Figure 3. Wc present the evolution of tlie scalar field amplitude for different values of the scalar field 
mass measured by a detector located at f 000 in units of the mass of the black hole, for the / = 
case. We found that the decay rate increases with the mass of the scalar field. 



In Fig. 4 we show the time dependence of the density of the scalar field p, in fact the 
physical density, that is, constructed with the physical and not the conformal scalar field. We 
again assume the density has a behavior of the type p ~ ^/f, in order to obtain estimates 
on a modest polynomial-like behavior in time. Other possibility we consider here, is the 
non-trivial self-interacting parameter A. In Table 1, we present the tail decay exponent for 
different values of the scalar field mass and A. In this table we also show the massless 
case for comparison, and we can notice that the effects on the polynomial exponent is more 
important when the mass is changed than when the self-interaction term is included. We 
consider values of A = ±1.5. 



A \ m| 


0.0 


0.1 


0.2 


-1.5 


p = -2.1179 


p = -4.0478 


p = -6.0397 


0.0 


p = -2.2412 


p = -4.0478 


p= -6.0398 


1.5 


p = -2.3952 


p = -4.0481 


p = -6.0389 



Table 1. In this table, we show the tail decay exponent for different values of the scalar field mass 
and the self-interacting parameter A measured by a detector located at 1000 for / = 0. As we can 
see, the effects that A introduces are minimal for the window of parameter here studied. In all the 
cases our fittings have an uncertainty smaller than 0.2%. 



Finally we show in Fig. 5 the convergence and self-convergence of one of our simulations, 
showing that our algorithms work properly. The order of convergence Q is consistent with 
our numerical implementation of sixth order in space {Q = 6) and fourth order (Q = 4) of 
the time integrator of the method of lines. Notice that the resolution in time is increased in 
the same proportion as the resolution in space, thus showing the convergence of the time and 
spatial implementation. The convergence factor appears closer to 6 than to 4 because the 
spatial approximation dominates over the time integration error, which is explained because 
after the profile of the scalar field has fiattened we are evolving nearly constant functions. 

Astrophysical scenarios and other models. In order to apply our numerical findings 
to astrophysical scenarios, we consider a supermassive black hole with mass IO^Mq, and 
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Figure 4. On the left we show the behavior of the physical energy density of the scalar field p 
normalized with the maximum of the density at initial time, as measured by detectors located at 
various distances from the horizon for I = 0,1,2 and m?g = 0.1. On the right we show a polynomial- 
like exponent that assumes the density behaves as ^ X/f. We use these bounds in order to estimate 
the dilution rate of the scalar field energy density in astrophysical situations. 



calculate the time it takes to the energy density of the scalar field to decrease five orders 
of magnitude for various values of I. We summarize our estimates in Table 2, where for 
completeness we have included the cases tub = and = 0.1. We again construct this 
table considering the polynomial decay stage only. 

The results indicate that the density decreases five orders of magnitude approximately 
every tenth of a year, independently of the initial amplitude of the scalar field, it would 
therefore dilute very quickly and would not survive cosmic time scales. This result is strongly 
related to the spherical symmetry that has been assumed, because the cross section for the 
accretion is maximal. 

Nevertheless, more realistic models assuming spherical symmetry that consider the dark 
matter to be an ideal gas show different scenarios, for instance, when the gas is pressure-less 
it was found that within cosmic time-scales the mass of the black hole shows an exponential 
growth, whereas when a non-zero pressure is assumed, the pressure itself contributes to 
stabilize the accretion into the black hole and provides accretion mass rates more consistent 




Figure 5. We show the convergence tests for one of our runs, in this case for ~ 0.1 and 
I ~ and resolutions Aa;i = 7.6 x 10~^, Ax2 = Axi/2, Axs = Axi/4. In the left panel we 
show the order of convergence Q of the L2 norm of the constraint C = — dr4> to zero, defined by 
L2{C)[using AX2] / L2{C)[using Ax^] = 2'^ because the ratio between resolutions is 2. In the right 
panel we show the order of self-convcrgcnce Qsc of 0, using the L2 norm of the differences between 
the value of the scalar field for the various resolutions. 

with observations [17]. In the present paper however, we do not measure accretion mass rates 
and actually focus on the considerable decrease of the scalar field density in time, because 
the scalar field propagates also toward scri +, whereas in the ideal gas case the particles of 
the fluid show speeds smaller than that of light. 

Another property of the scalar field as dark matter is that a phase space analysis seems 
to be less practical than the solution of the Klein-Gordon equation itself; such phase space 
analysis would be helpful to study an adequate cross section with the black hole, and provide 
more realistic conditions for the scalar field to propagate, like rotation and non-zero impact 
parameter, compared to the spherical symmetry we present here, and also more similar to 
standard studies of accretion of dark matter into supermassive black holes [15]. 

5 Conclusions 

The rapid decay of the energy density of the scalar field for the case of a super-massive 
black hole, indicates that the scalar field cannot be maintained around the black hole dur- 
ing cosmological time scales in the whole space; either the scalar field is being accreted or 
escapes through future null infinity. The fact is that when a Schwarzschild black hole -that 
is asymptotically flat- is considered to be a black hole candidate, scalar fields tend to vanish 
from the spatial domain. 

One potential ingredient that would help at maintaining massive scalar field densities 
during longer times is the rotation of the black hole and also of the scalar field. It would 
be of the major interest to use foliations that approach future null infinity and study if the 
same effects occur and also the study using scalar field configurations with non-zero angular 
momentum. In fact our results assume the maximum cross section of accretion due to the 
spherical symmetry we assume, which in turn works as upper bounds accretion rates in more 
general cases. 

Yet another possibility is to consider solutions that asymptotically containing cosmolog- 
ical constant, which would be appropriate if a background energy density in the universe is 
assumed. This would imply that black hole candidates should not be considered to be asymp- 
totically flat. Other possibilities may include black hole candidates of a different nature like 
boson stars [25]. 
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1 X 10-25 - 1 X 10-30 


0.137 
0.075 


m| = 0.1 


I 


f [M] 


At [M] 


pi-pf 


At [years] 





1000 

58 


905 
628 


1 X 10-2i - 1 X 10-2« 

1 X 10-24 - 1 X 10-29 


0.14 
0.09 


1 


1000 

58 


680 
442 


1 X 10-2^ - 1 X 10-2'^ 

1 X 10-25 - 1 X 10-30 


0.10 
0.06 


2 


1000 

58 


567 
319 


1 X 10-25 - 1 X 10-3U 
1 X 10-28 _ 1 X 10-33 


0.08 
0.05 



Table 2. We consider the Schwarzschild solution to have the mass of a supermassive black hole with 
mass M = 10^ Mq. We calculate the time in years that takes the physical p to decay its value five 
orders of magnitude, which we consider to be enough to decide on the possibility that the scalar field 
could remain around the black hole for time scales of cosmological order (gigayears). We consider 
three cases / = 0, 1,2, and measure the decay time using two detectors located at different distances 
from the horizon {f[M]). We locate a time interval At[M] that takes the density to decrease an 
amount pi — pf and calculate it in years At[years]. 
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A Appendix 

In this appendix we show the construction of equations (3.3) starting from (3.1). We start 
with the D'Alembertian operator for a general metric written as 



In our spherically symmetric metric, ^/—g = a^r"^ sin 6 is the determinant of the metric. 
Then (3.1) reads 




(A.l) 




1 



R<t)T - {mn~'^cl)T + A0|) , 



(A.2) 



6 
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which imphes 



T - -R(I)T - {mil. "^(pT + >^(Pt) 



d. 



For our first order variables 



+ 



+ 



1 

1 

1 



A- 



Xdt4)T + ^drCpT 

a a 



sin ( 



deismOdecpT) + -T—^d^{d^(t)T 



sm ( 



1 



R(j)T - {mQ,~'^(j)T + X4>t) ■ 



a a 



the above equation takes the form 



(A.3) 



(A.4) 
(A.5) 



dtTT = 



sm 



sm ( 



cP{t,r) 



Rcl)T + (mO + A(/)|) 



(A.6) 



where we have assumed 4)T{t.,r,9,ip) = 4>{t,r)Yim{9,ip). Now, the expression between brack- 
ets in the second term is the definition of the Laplacian of spherical harmonics, which can 
be expresed by the following relation 



sm 



sin( 



-lil + l)Ylmi9,ip). 



(A.7) 



Finally, we obtain the evolution equation for vr and the other two evolution equations for ^/^ 
and (pT can be obtained easily from the expressions (A.5) 



dtlp = dr ^^TT + Pip 



OL 

dt4>T = -vr + p-ip. 

7 



(A.8) 
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